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Simulating leaf net C0 2 assimilation rate of C 3 & C 4 plants and its re¬ 
sponse to environmental factors 
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Abstract: Basic structure and algorithm of leaf mechanism photosynthesis model were described in first part of 

this study based on former researcher results. Then, considering some environmental factors influencing on leaf 
photosynthesis, three numerical sensitivity experiments were carried out. We simulated the single leaf net C0 2 
assimilation, which acts as a function of different light, carbon dioxide and temperature conditions. The relation¬ 
ships between leaf net photosynthetic rate of C 3 and C 4 plant with CO 2 concentration intercellular, leaf tempera¬ 
ture, and photosynthetic active radiation (PAR) were presented, respectively. The results show the numerical ex¬ 
periment may indicate the main characteristic of plant photosynthesis in C 3 and C 4 plant, and further can be used 
to integrate with the regional climate model and act as land surface process scheme, and better understand the 
interaction between vegetation and atmosphere. 
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Introduction 1 

It is well recognized that the global climate is be¬ 
coming warmer and warmer owing to increased 
emission of greenhouse gases affected by human 
activities. The better method of predicting future cli¬ 
mate trend is to use regional and global climate 
model. Recent years, a series of such kinds models 
have been developed to describe the interaction be¬ 
tween the biosphere to atmosphere, and several of 
land surface interaction complex process schemes 
have been incorporated into regional or global cli¬ 
mate models (Dickinson et at., 1986; Sellers et al., 
1986, 1992, 1996; Xue et al. 1991; Foley et al, 1996). 
As we know, the land surface process in regional 
climatic model (RCM) employed in eco-physiology 
are designed to describe the basic photosynthesis 
process at the leaf, canopy and regional levels to¬ 
gether with the stomatal conductance, vegetation 
transpiration, C0 2 flux change within leaf. Repre¬ 
sentation of C0 2 exchanges at the sur¬ 
face-atmosphere interface is an important challenge 
for assessing the impact of climatic change on the 
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surface energy and water budget, and it has been 
suggested that the exchange of water vapor and C0 2 
between vegetation canopies and the atmosphere is 
strongly controlled by the physiological processes 
governing photosynthesis and stomata conductance. 
Thus, developing models to estimate photosynthesis 
is certainty. 

The purpose of current study is to simulate the 
response of photosynthesis of C 3 & C 4 pathway 
plants to environmental affecting elements based on 
former researcher schemes; additionally, is to under¬ 
stand the realistic photosynthesis process of differ¬ 
ence photosynthesis pathway in vegetation. Fur¬ 
thermore, to develop a more realistic mechanize 
photosynthesis scheme for integrating RCM. 

Three pathways of plant photosynthesis 

Generally, based on the photosynthesis pathway, 
the natural flora can be divided into three categories 
(Monson et al. 1989), namely C 3 plant (by all trees 
and many herbs, C0 2 assimilation through reductive 
pentose phosphate cycle), C 4 plant (tropical herbs 
and warm grasses; 20 family; >12 000 species, C0 2 
assimilation through PEP carboxylase), and CAM 
(crassulacean acid metabolism) plant (carnification 
plant; 26 family; >500 species). Tablel lists the main 
photosynthesis characteristics of C 3 , C 4 and CAM 
plants. In present study, numerical algorithm will be 
adopted to approach some important photosynthetic 
characteristics of C 3 and C 4 plants. 
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Table 1. C 3 , C 4 and CAM plants photosynthesis and other eco-physiological characteristics 


Characteristics 

C 3 plant 

C 4 plant 

CAM 

First step, Catalytic Reaction 

Rubisco E 

Rubisco E 

Rubisco E 

Fixed C0 2 

C0 2 immobilization enzyme 

RuBP E 

PEP carbonlase 

PEP E 

Representative plant 

Wheat (crop), Pine (tree) 

RuBP E 

Corn (crop), Millet (cereal) 

RuBP E 

Cactus (shrub) 

Leaf structure 

Bundle-sheath (undeveloped 

Bundle-sheath (developed) 


Rmax (mgC0 2 - dm' 2 - h' 1 ) 

Low (15-40) 

High (35-80) 

Low (1-4) 

C0 2 compensation point (mg- kg' 1 ) 

High (40-70) 

Low (0-10) 

Dark period (0-5) 

Light respiration 

High 

Low 

24-h (0-200) 

Low 

Light saturation point 

Low ( =1/4-1/2 Rmax) 

High (above R max ) 

Unfixed 

Suitable temperature of Photosynthesis 

Low (15-25) 

High (30-47) 

-35 

Suitable temperature of growing 

Low 

High 

Widely 

Resistance to arid climate 

Weak 

Strong 

Changed 

NPP (g- m' 2 - d' 1 ) 

Low (19.5+3.9) 

High (30.3±13.8) 

Changed 

NPP (t- hm' 2 - a' 1 ) 

Low (22.0±3.3) 

High (38.6116.9) 

Changed 

Water demand (g water/g dry-matter) 

Large (450-950) 

Small (250-350) 

Very small (50-55) 

Increasing of dry-matter /Doubling C0 2 

high (>20%) 

Low (<10%) 

Unfixed 


The leaf photosynthesis scheme 


and for C 4 plant 

w c ~ U m 


(3) 


Modeling method 

After 1980s, some mechanism photosynthesis 
simulating models have been set up and developed. 
The significant progress was Farquhar (1980) and 
Collatz (1990,1991) works. In their models, the pho¬ 
tosynthetic rate (A) is a function of PAR, leaf tem¬ 
perature, C02 flux within the leaf, and the Rubisco 
enzyme capacity for photosynthesis (Sellers et al. 
1996a). In both cases the rate of gross leaf photo¬ 
synthesis, A, is calculated in terms of three potentially 
limiting factors. In this paper, the term A of C 3 & C 4 
plants are similarly modeled as the minimum of three 
potential capacities, be expresses as 

A < min (w c , w e ,w s ) (1) 

where, A is gross Leaf photosynthetic rate 
(mol-m' 2 s'); w c is represents the rate of gross pho¬ 
tosynthesis when the leaf photosynthetic enzyme 
(Rubisco of leaf enzyme) was limited; w e is light- lim¬ 
ited rate of gross assimilated rate (mol-m^-s' 1 ); w s is 
the limitation associated with export of the photo¬ 
synthetic products (for C 3 vegetation), or 
PEP-Carboxylase (for C 4 plant) limitation on photo¬ 
synthesis (mol-m' 2 -s' 1 ). 

For the physiological limit of photosynthetic en¬ 
zyme on assimilation rate, u c can be written as: 
for C 3 plant 




c, -r* 


C,+^ c (l + 0 2 /^o)J 


where, w c is Rubisco-limited rate of assimilation 
(mol-m'^s' 1 ); l/ m refers to maximum rate of catalytic 
capacity of Rubisco (mol-m'V 1 ) at T c ; Q is partial 
pressure of C0 2 in leaf interior(Pa); 0 2 is partial 
pressure of 0 2 in leaf interior (Pa); F* is C0 2 com¬ 
pensation point (Pa) = 0.5O 2 /S; S is Rubisco speci¬ 
ficity for C0 2 relative to 0 2 (=2600‘0.57 a ); Kc is 
Michaelis-Mnten constant for C0 2 (Pa), and its value 
is equal to 30x2.1 01 ; is inhibition constant for 
0 2 (Pa), and its value is 30 000x1 .2 a , where Q, is Q 10 
temperature coefficient, equal to O.lx (T c -25), T c is 
the leaf temperature (°C). 

For C 3 plant, the term V m is given by 

V = T -_ (4) 

{1 + exp[0.3(T c — 38)] } 

and for C 4 plant, it is given by 

y __ ^ma \Qt _ /g\ 

m ” ^+ exp[0.3(T c -40)]}^ + exp[0.2(5-r c )} 
where, Umax is maximum leaf catalytic capacity of 
Rubisco at 25°C, the values V ma x were given by 
1x10' 4 (mol-m'? s' 1 ) for C 3 broadleaf deciduous tree 
and 3x1 O' 5 (mol-m' 2 s' 1 ) for C 4 grass (Sellers et al. 
1996b). 

The light-limited rate of assimilation w e is given by 
(Foley etal. 1996): 

w e = 0.068/ PAR [ Cl * ■] for C 3 plant (6) 
C; +2L* 


( 2 ) 
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w e =0.0332/,^ for C 4 plant (7) 

where, l par is the incident PAR on the leaf 
(mol m' 2 s _1 ), and its value depends on the assumed 
quantum efficiency for C0 2 uptake and the leaf scat¬ 
tering coefficient for PAR. In this study, it can be 
written the PAR multiplies an absorbed factor by 0.86 
(Collate ef a/. 1991). 

A third limiting rate has been defined for C 3 and C 4 
photosynthesis by Collatz et at. (1991, 1992), re¬ 
spectively, w s is viewed as the capacity for the export 
or utilization of the products of photo synthesis in the 
case of C 3 and as the C0 2 -Iimited capacity for C 4 
photosynthesis. 

w s = 0.5 V m for C 3 plant (8) 

w s = 2x1 0 4 V m c/p* for C 4 plant (9) 


response of C 3 & C 4 plants to PAR, leaf temperature, 
C0 2 flux within leaf for has been addressed. Fig. 1 
shows the dependence of A n in leaf level on PAR, 
reflects that the C 4 light saturation points is obviously 
higher than that of C 3 plant. 



0 500 1000 150 

PAR( Photosynthetic Active Radiation)(Wm-2) 


where, p* is atmospheric pressure (Pa). 

The actual rate of gross photosynthesis, A, is cal¬ 
culated as the smoothed minimum of above three 
limiting rates: 


P\w 2 p - w p (w c +w e ) + w e w c =0 
P 2 P 2 - P(w p +w s ) + w p w s = 0 


( 10 ) 


where, is coupling coefficients, w p is smoothed 
minimum of w c and w e (mol m^ s 1 ). 

The net assimilation A n , is then given by 

A n =A-R d (11) 

where, fl d is leaf respiration rate(mol-m’ 2 s' 1 ), Collatz 
et al. (1991, 1992) scaled fl d to the leaf carboxylase 
Content by 

*d=/# m (12) 


where, f d = 0.015 for C 3 plant, and f a = 0.025 for C 4 
plant. 


Fig.l The relationship between leaf A n and PAR for C 3 
and C 4 plants 

Fig. 2 presents the relationship between A n and 
leaf temperature. In X-axis, before 20°C, with leaf 
temperature increases, the increasing trend of C 3 
plant A n is higher than that of C 4 plant, however, after 
20°C, there is a significant increasing of A n for C 4 
plant. 
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The numerical experiments and some results 

In this study, two plant types were selected as ex¬ 
perimental plant species, the temperate deciduous 
broadleaf tree (C 3 plant), located in north part of 
Changjiang River, and maize (C 4 plant), which is 
grown in North China Plain. Some experiments test¬ 
ing were performed via considering the some factors 
influencing on leaf photosynthesis. Here, we consider 
three kind of effect situation on photosynthetic rate. 
One is the dependence of leaf photosynthesis on 
Photosynthesis Active Radiation (PAR), second is the 
relationship between leaf photosynthesis and C0 2 
concentration intercellular, third is the relationship 
between leaf photosynthesis and leaf temperature. 

We simulated the net C0 2 Assimilation (A n ) in leaf 
level, the model at the leaf level can also be used to 
modeling A n in the larger scales by scaling up. A n 


Fig.2 The relationship between leaf A n and leaf 
temperature for C 3 and C 4 plants 

Anyway the maximum values of A n (C 4 plant) 
nearly against 36°C. But in case of C 3 plant, the leaf 
temperature just less than 30°C when A n increases 
its peak. The result implies a better elucidation for C 4 
plant originated from tropic zone, to some extant, also 
reflects the light saturation point of C 3 plant is less 
than that of C 4 plant (Harold et at. 1978). 

Fig. 3 & Fig 4 show the relationship between A n 
and Cj about C 3 and C 4 plants. Owing to C0 2 com¬ 
pensation point of C 3 plant is higher than that of C 4 
plant, so the Fig. 3 presents there is a linear positive 
correlation between A n and Q. Nerveless, in Fig.4, A n 
of C 4 plant is close to saturation with the Q increasing. 
This result also reveals the increasing of dry-matter of 
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C 3 plant is higher than that of C 4 plant under C0 2 

doubling. References 
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Fig. 3 The relationship between net leaf A n and Ci for C 3 
plant 



Fig. 4 The relationship between leaf A„ and Ci for C 4 
plant 


4. Conclusion 

Preceding the numerical results present the de¬ 
pendence of leaf photosynthesis rate (An) on Photo¬ 
synthetic Active Radiation (PAR) curve, the leaf net 
photosynthesis ( A n ) vs. leaf temperature (T c ) re¬ 
sponse curves in C 3 and C 4 plant, and the leaf net 
photosynthesis (A n ) vs. intercellar C0 2 concentration 
(Ci) response curves in C 3 and C 4 plants. The model 
presented here provides a plausible mechanism for 
feedbacks of photosynthesis to multiple environ¬ 
mental forcing such as PAR, T c and Q. Simulation 
indicates that the environmental factors have signifi¬ 
cant influences on the leaf net photosynthetic rate. 
Whereas, the additional observed experimental 
studies at the leaf level should be required to test the 
mechanism proposed here in future work. 

To arise from this study, we have demonstrated 
that the numerical experiments may reveal the main 
characteristic of plant photosynthesis in C 3 and C 4 
plant, and further seem to be used to integrate with 
the regional climate model and act as land surface 
process scheme, and better understand the 
interaction between vegetation and atmosphere. 
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